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Abstract 

We study numerically the gravitational field of a star made of massive and neutral 
string states for the case in which the dilaton is massive. The solution exhibits very 
simple scaling properties in the dilaton mass. There is no horizon and the singularity 
is surrounded by a halo (the physical size of which is inversely proportional to the 
dilaton mass) where the scalar curvature is very large and proportional to the square 
of the dilaton mass. 



1. Neutral and scalar string excitations which are massive already in ten di- 
mensions [m represent the simplest form of that exotic state of matter in which 
the constituents are string states not belonging to the massless sector of the string 
spectrum to which instead ordinary matter, like quarks and leptons, presumably 
belongs (their masses originating through some still unknown symmetry breaking 
mechanism) . 

The physical relevance of such string matter stems from several studies which 
indicate that, at high energy density, the most probable configuration in string 
theory is the one in which most of the massive states are excited. According to this 
view, an electrically neutral collapsing star of sufficiently large mass would start 
out as a celestial body made of ordinary matter but would eventually evolve into a 
string star in which most of its mass is now carried by neutral string excitations. 

It is therefore of some interest, we believe, to study the gravitational field around 
such a string star. 

The scalar, neutral and massive string excitations couple only to gravitons and 
dilatons; therefore, the gravitational field around a star made of a large number 
of such states is described by Einstein's equations in which the energy-momentum 
tensor is the one for a scalar field. In the original string scenario, these equations 
arise as conditions for the vanishing of the one-loop beta functional required by 
Weyl invariance; the dilaton is massless and the field equations admit an exact 
solution [|, ^ . This solution, and its relationship to string theory, has been discussed 
recently in ||^, to which this letter is closely related. The most significant feature of 
such a solution is the absence of horizon. The relevant elements Qrr and Qqq of the 
space-time metric are shown in fig.l, where they are compared to the Schwarzschild's 
ones. As it is possible to see from the figure, while the two solutions are equivalent at 
large distances, they are remarkably different closer to the gravitational radius: at 
the horizon Schwarzschild's grr diverges and —goo, which is just its inverse, crosses 
into negative values, whereas the ones corresponding to the string star solution 
are never either negative or infinite. These results seem to be in agreement with 
previous work in which the existence of a static and uncharged static solution 
with a scalar field has been shown to be incompatible with the presence of a non- 
singular horizon (see, also, 0). 

It is tempting to speculate that also in a realistic scenario, derived from string 
theory, in which ordinary matter, after obtaining a mass much lower than Planck 
mass by some yet unknown mechanism, continues to behave with respect to the 
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dilaton as it were in the massless sector of the string spectrum. According to this 
view, the dilaton field remains massless as well as decoupled to ordinary matter at 
large distance. Such a possibility can be entertained without violating any of the 
present experimental observations. In this case all the results of would hold true. 

In this letter we address the question of what happens to the solution with a 
massless dilaton in the opposite case in which the dilaton does couple to ordinary 
matter and therefore — because of the current astronomical bounds [0 on Newton's 
inverse- square law — it has to become massive. Such an investigation is, for its 
very nature only preliminary, because of our present ignorance of how the basic 
ingredients of the superstring, like world-sheet conformal invariance, determine the 
structure of the effective low-energy theory with a massive dilaton (see for a 
related discussion in cosmology). 

The reader, confronted by a growing literature on the subject of solutions to 
graviton plus dilaton (and, plus Maxwell fields) gravity in four space-time dimen- 
sions, should bear in mind that the general solution for a source with an arbitrary 
coupling to massless dilaton, graviton and Maxwell fields is contained and discussed 
in the papers of ref. 0. Paper [^] discusses the special case in which the coupling 
to the massless dilaton field is a function of the Maxwell charge and it vanishes as 
the charge goes to zero. Our solution is the one in which the source is neutral with 
respect to the the Maxwell field but nonetheless has a non-vanishing coupling to 
the dilaton field. 



2. The field equations are easily written as variations of the action 
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We have restricted ourselves to the static and spherically symmetric case and there- 
fore used the parametrization 
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for the space-time metric, the signature of which is — H — h +; prime means differ- 
entiation with respect to the radial coordinates, (p = K,(f is the rescaled dilaton field 
ip and = StiGn, Gn being Newton's constant. 
In general, 

W{^) =(f)^ + a(f)^ + p^^ + ... , (4) 

where the first term gives rise to the mass, and the remaining ones to an effective 
potential for the dilaton field. 
Two questions arise. 

The action (|T]) without the mass term is obtained from the usual one 

S = ^ f dxv^e-2* [7^ + 4V0 ■ V0] (5) 

by the rescaling g e^^g and -\/20/2, and it represents the low-energy effective 
theory one derives from string theory. After symmetry breaking, other terms beside 
the mass m^0^ may appear and give rise to a potential for the dilaton field as in (^. 
Consider, for instance, a term like Even though the coupling to ip^ is typically 
of the order of Gn^"^ (that is, (3 of order of one), the net contribution to the 
equations of motion is of the same order as the ones included in (|l]), as one can 
readily see by inspection. 

We have explicitly checked our solution against such terms in the potential and 
verified that its nature, qualitative features as well as most of the numerical values 
are left unchanged for values of the parameters a, (3 up to order of ten. The 
only case in which there is a substantial change in behavior is when we take, for 
example, a > and j3 and the coefficients of the other even power terms vanishing, 
thus forcing the potential to become unstable (recall that is negative). Therefore, 
except for this somewhat pathological case, we are entitled to disregard for the sake 
of simplicity these extra terms, and take in what follows W{(j)) = </)^ only. 

A second question concerns possible string corrections to (|l|). These are, for 
example, of the type a'TZ^, where a' is the inverse of the string tension. However, 
in this case, one finds that even for very large values of the scalar curvature (as we 
are in fact going to have) such terms are suppressed by a factor a'm^ with respect 
to the term linear in the curvature TZ, and can accordingly be neglected altogether. 

3. Equations @ do not seem amenable of any exact solution. We study them 
by numerical methods. 
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It is convenient to recast (0) as a system of four first-order differential equations: 

1 y2yl 
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y'z = ^^p^vm 

y'. = (6) 

by means of four dimensionless functions: 

Z/i = 

ys = pe(^-^)/2 

y. = p'e^'^^'^'^' . (7) 

In (D), p = r/rg, where Vg = 2GnM is the gravitational radius of the star. The 
prime is now differentiation with respect to the coordinate p. 

The parameter ^ = mrg is very large for the physically interesting scenario 
in which the Compton wavelength of the dilaton 1/m is much shorter than the 
gravitational radius; we have representatively considered values of C, between 10 to 
10^0. 

We set the boundary conditions at large distance for (||) as follows: 
Z/2 1 

1/4 p-1, (8) 

that is, we require that both goo and grr must go into Schwarzschild metric for 
sufficiently large distances. 

Less clear is the asymptotic value for the dilaton field (p. The exponential decay 

l/i-^0oe-«7P, (9) 

that we expect by solving the corresponding equation, leaves the constant (^o un- 
determined. While for the massless case this indeterminacy can be resolved by 
comparison with string perturbation theory P], the problem seems more compli- 
cated here. Luckily, as we shall see, the solution we find is very stable over a large 
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range of changes in the initial condition for the dilaton field. We therefore fix it by 
assuming in the region p — 1 <^ 1 that 



y^^Ko(2i^p-lJ , (10) 

as it is suggested by the study of the dilaton equation in a fixed background, and 
putting by hand the numerical coefficient in front of the Bessel function. This choice 
makes the dilaton field behave as in (|^) for p — 1 ^ 1, and as 



exp -2^Vp - 1 

^ — oc I — r — 



and 



yi^log 2Wp-1 (12) 



for, respectively, ^ p — 1 ^ 1 and < p — 1 ^ 

This procedure amounts to placing the position at which the dilaton field begins 
to be different from zero. Within an acceptable range of initial values, once it has 
been fixed, the dilaton field rapidly goes into its asymptotic (and universal) behavior 
beyond the gravitational radius. 



4. We can thus solve the system of first order differential equations (|^) as an 
initial- value problem by means of a standard program based on a variable-order. 



variable-step method implementing the Backward Differential Formulas. 

The results of the numerical integration are displayed in figs. 2a-2b for the 
metric and 3 for the dilaton field, where they are plotted for the special value ^ = 10. 
Although such a small value is not physically relevant, it is useful because it makes 
the plotting more readable. Moreover, the solution exhibits very simple scaling 
properties that allow to infer from the plot of figs. 2a-2b and 3 the corresponding 
ones for higher, and more realistic, values of ^. 

As we have already pointed out, we obtain a solution which is surprisingly stable 
with respect to variations in the dilaton mass, the asymptotic value of and the 
parameters in 1^(0). The features of the solution we would like to emphasize are 
that: 

• like in the massless case, the horizon is absent; 
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• there is a very narrow region near the would-be horizon where there is a very 
sharp transition away from Schwarzschild's solution; 

• outside this region — that is, for r > Vg — the field can be neglected and the 
solution is indistinguishable from the Schwarzschild's one; 

• inside this region — that is, for r < Vg — the solution becomes stable and de- 
pends very little on details. 

By varying the input value of C, and comparing the outputs so obtained we have 
verified that the behavior of the space-time metric for p ^ is the following: 

grr ^ ^POP^O/^^^ 

-^00 ^ ^5op^o/2-i (13) 
while the dilaton field remains unchanged and: 

0^gologp-Co. (14) 
Also the scalar curvature behaves in a simple manner: 

7^^e'|^p-^'+'^°/'^ + 0(logp)| . (15) 

The numerical values for Qo, Co and Pq seem to be universal, independent of the 
dilaton mass and every other parameter; they turn out to be Qq ~ 2.5, Cq ~ —1 
and Pq ~ 4. 5*0 seems to depend more on the input conditions, however being 
always of order one. Notice that the solution for the massless case was of the 
same kind with Ql = 2(^2 + 1)7(^2 - 1). 

The behaviors ([T3|)-(p!5D are well satisfied inside the Schwarzschild radius, whereas 
away from it the solution is more sensitive to the initial conditions and does not 
follow a simple scaling law in ^. 



5. Let us briefly discuss the solution we have found. It resembles in a quali- 
tative manner the massless case. Away from the would-be horizon, the metric is 
indistinguishable from Schwarzschild's. However, when distances of the order of 
1/.^^ are reached, the solution departs dramatically from Schwarzschild's and the 
horizon disappears. For large values of ^, this change takes place in a very thin shell 
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around r^. This can be compared with the massless case in which the departure 
occurs at distances of the order of Vg instead. 

At the same time, the scalar curvature (see fig.4) grows in a very steep manner 
and produces a high- curvature sphere around the singularity. This is to be expected 
since in this region (j) changes from being vanishingly small to being of order one, 
thus acting back on the space-time metric and giving a large contribution to the 
curvature as well. Fig.4 makes clear how the massive case differs from the massless 
one only in that the growth of the curvature is more abrupt (again, it takes place 
in the region of size 1/,^ before r^) and steeper. They become equivalent closer to 
the singularity where the dilaton field becomes effectively massless in any case. 

The physical size of the halo of large curvature is 



and therefore of order 1/m in standard units. 

We feel that these results arc very natural. The presence of a mass for the 
dilaton field, by making available a second parameter to the problem, changes the 
scale at which the departure from Schwarzschild solution takes place and also the 
magnitude of such departure; it does not however change the nature of it, that 
remains the one found in the massless case. 
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Captions 



Fig.l: Schwarzschild's space-time metric compared to the string star's for a massless 
dilaton field. 



Fig. 2a: Schwarzschild's goo compared to the string star's for a massive dilaton field. 



Fig. 2b: Schwarzschild's Qrr compared to the string star's for a massive dilaton field. 



Fig. 3: The massive dilaton field. 



Fig.4: The scalar curvature: Massless vs. massive dilaton field. 
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